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Abstract 

We analyze the effect of an external electric field on the electronic structure of molecules which 
have been recently studied as molecular wires or diodes. We use a self-consistent tight binding 
echnique which provides results in good agreement with ab initio calculations and which may be 
applied to a large number of molecules. The voltage dependence of the molecular levels is mainly 
linear with slopes intimately related to the electronic structure of the molecules. We emphasize 
that the response to the applied voltage is an important feature which governs the behavior of a 
molecular device. 



I. INTRODUCTION 



Recent measurements of single molecule transport properties^ - - represent important con- 
tributions to the eventual realization of molecular electronics. This is accompanied by an 
increasing theoretical effort to understand the relationship between the electronic structure 
of the molecules and the current-voltage /(V) characteristics. Until recently, most of the 
calculations based on tight binding (TB) or ab initio methods^ did not explicitly include 
the effect of the applied voltage on the electronic structure of the molecules. This is a se- 
vere drawback since several volts can be applied to the electrodes in contact with a single 
molecule. In the some cases^, the energy of the molecular electronic levels is approximated 
by a linear function of the voltage but its slope is not calculated or is approximated by 1/2. 
An approach to overcome these limitations has been proposed^ recently on the basis of 
a Hubbard Hamiltonian^ and has been applied to the azulene molecule^ and to polyacene 
wiresS. But the Hubbard Hamiltonian is a model Hamiltonian which only gives qualitative 
predictions. On the other hand, ab initio methods cannot be easily applied to complex 
systems. Thus there is a need for simpler methods which correctly predict the response of 
the molecules and the dependence of the electronic levels as function of the applied bias. 
This is a prerequisite to calculate the I(V) characteristics of a molecular device. In this 
paper, we present self-consistent TB calculations of various molecules. We study the evolu- 
tion of the dipole and of the electronic levels with the voltage. We obtain results in good 
agreement with ab initio calculations in the local density approximation (LDA). This is a 
good test of our method and allows to extend it to large systems. We conclude that our 
self-consistent TB technique is a promising approach to study the transport properties in a 
molecular device and that it may be applied to a wide range of systems. 



II. CALCULATIONS 



The molecules considered in this work are shown in Figure [TJ In order to compare 
with previous results^-, we study tetracene and azulene molecules as prototypes of molec- 
ular wires. The calculations are also applied to quinolinium tricyanoquinodimethanide 
(Q-3CNQ) and to 3,5-dinitrobenzyl 7-(l-oxohexylamino)-pyren-2-ylcarbamate (OHAPy-C- 
DNB) which are respectively D-7T-A and D-cr-A molecules. D and A are respectively elec- 
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tron donor and acceptor. 7r and a are respectively "pi" and "sigma" bridges. These two 
molecules are intensively studied to make molecular diodes^. The electronic structure of 
all the molecules is obtained in LDA and in TB. In LDA, we use the DMOL code— with a 
double numerical basis set (two atomic orbitals for each occupied orbital in the free atom) to- 
gether with polarization functions (2p for H, 3d for N, O and C). The exchange-correlation 
energy is approximated by the density functional of ref^. The self-consistent TB tech- 
nique is presented in another publication 13 - where we calculate the electronic structure of 
thienylenevinylene oligomers, showing a good agreement with LDA calculations and exper- 
iments. C, N, O atoms are represented by one s and three p atomic orbitals, and H atoms 
by one s orbital. The non-diagonal terms of the Hamiltonian matrix H are restricted to first 
nearest-neighbor interactions and to two-center integrals. They depend on the interatomic 
distance following Harrison's rules^. The self-consistency is incorporated in the diagonal 
terms H ia ^ a (a is the orbital index, % is the atomic index) 

H iaiia = H® aia - i ^ J = - eV ext (i) (1) 

where Hf a ia define the atomic levels, U is the intra-atomic Coulomb energy, is the 
distance between atoms i and j, and Q., is the net charge on the atom j. The TB parameters 
are determined by a fit of the LDA electronic structure of simple molecules^ 3 -. The parameters 
for C-C and C-H interactions are given in Ref^ 3 -, those for C-N, N-H, C-O, N-0 interactions 
are in Table [H The selfconsitency is obtained with an usual iterative method which is 
converged when the atomic charges Qj in eq. ([T]) correspond to the charges calculated from 
the eigenstates of the Hamiltonian. V ex t in eq. ([[]) is the electrostatic potential resulting from 
the applied electrostatic field which, for simplicity, is assumed homogeneous, corresponding 
to a situation where the electrodes would be far from the molecule. Obviously, to calculate 
the I(V) curve, we should consider the chemical and electrostatic interactions between the 
molecule and the electrodes. This can be done in TB but this is beyond the scope of the 
present paper. 
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III. RESULTS AND DISCUSSIONS 



The applied electric field E is parallel to the long axis (z) of the molecules (Figured]). 
The origin of the electrostatic potential (V e xt = ) is defined at 2 A from the left side 
of the molecules. Thus V ext (i) = Ezi where z^ is the coordinate of the atom % along the 
axis z. All the results presented in this paper are plotted as a function of the electrostatic 
potential at 2 A from the right side of the molecules. Equivalently, a bias of 1 V corresponds 
to an applied electric field of 88 mV/ A, 64 mV/A, 57 mV/A and 46 mV/A respectively 
for azulene, tetracene, Q-3CNQ and OHAPy-C-DNB (the inter-electrode distance is 11.4 A, 
15.5 A, 17.4 A and 21.6 A, respectively). We plot in Figure 2 the voltage dependence of 
the dipole moment of the four molecules. Azulene has a small permanent dipole moment 
resulting from its nonalternant character—. Q-3CNQ is characterized by a large dipole 
moment due to charge transfer between acceptor and donor sites. OHAPy-C-DNB have 
weaker acceptor and donor characters. We see that the agreement with LDA is good. In 
particular, the overall magnitude and the slopes of the curves (or the polarizabilities) agree 
well. It means that the charge transfers induced by the applied bias are correctly described. 
In ref&, the dipole moment of azulene was not calculated, but equivalently the charge on 
the five-membered ring was plotted as a function of the voltage. Between -3 V and +3 V, a 
variation of this charge of 0.5e is reported to be compared with our much smaller value of 
0.3e. We present in Figure [3]the voltage dependence of the energy levels of tetracene. A very 
good agreement is obtained with LDA results, which was not the case with the Hubbard 
Hamiltonian^. The variation of the levels is mainly linear, with a slope 1/2 as expected from 
the symmetry of the molecule. In spite of an asymmetric atomic structure, quite similar 
results are obtained for Q-3CNQ (Figure HJ) with a slope 0.52 for the highest occupied orbital 
(HOMO) level and 0.49 for the lowest unoccupied orbital (LUMO) level in TB (resp. 0.56 
and 0.55 in LDA). This is a consequence of the derealization of these orbitals over the whole 
molecule, which results from an efficient coupling between acceptor and donor sites through 
the 7r bridge. The situation is completely different in the case of OHAPy-C-DNB (Figure 
[5]) because the a bridge leads to a much smaller coupling between the acceptor and donor 
parts. Therefore, the HOMO is mainly localized on the acceptor site and the LUMO on the 
donor site. As the donor site is closer to z = where the origin of the potential is defined, 
the voltage dependence is smaller for the HOMO than for the LUMO. All these important 
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tendancies are predicted similarly in LDA and TB, even if TB calculations underestimate 
the HOMO-LUMO gap compared to LDA in the case of OHAPy-C-DNB. Thus, even if Q- 
3CNQ and OHAPy-C-DNB are donor-acceptor molecules, their electrical characteristics will 
completely differ because of different voltage dependence of the levels. We have also studied 
the effect of an applied voltage on other molecules like thienylenevinylene oligomers^ and 
the 7-hexadecylquinolinium tricyanoquinodimethanide (Ci6H 33 Q-3CNQ, Figure [[]). The 
response is also linear between -3V and +3V, and the agreement between TB and LDA 
remains good. The case of Ci6H 33 Q-3CNQ is interesting because it differs from Q-3CNQ 
only by the replacement of an hydrogen atom by an hexadecyl group. Thus the HOMO and 
the LUMO are almost the same in the two molecules. However, their voltage dependence in 
CigH 33 Q-3CNQ is characterized by a much smaller slope (~ 0.2) due to the fact that a large 
part of the voltage drop takes place in the long Ci6H 33 group which is weakly polarizable, 
while the HOMO and the LUMO are localized in the remaining part of the molecule. Once 
again, we expect completely different /(V) curves for these two molecules in spite of their 
important similarities. 

IV. CONCLUSION 

We have presented a self-consistent TB method to study the response of molecules to an 
electric field. Our work justifies for a broad range of molecules the hypothesis of a linear 
dependence of the molecular levels with respect to an applied voltage^, even up to several 
volts. The self-consistent TB method provides a simple approach to calculate this depen- 
dence with enough accuracy. It represents a good compromise between model Hamiltonians 
and ab-initio calculations. It will be useful to simulate more complex systems, to calcu- 
late the charging energy of molecules^ and to predict the I(V) characteristics of molecular 
devices. 
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Interactions 


C-0 


N-H 


N-0 C-N 


O-H 


(ss a) 


-3.62 
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-3.32 -2.90 


-6.96 


(sp a) 


+4.96 


6.03 


+4.50 +3.68 
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(PP & ) 


+5.31 




+4.55 +5.23 




(pp 


-2.00 




-1.98 -2.00 




do 


1.54 


1.07 


1.54 


1.07 


Atomic levels 


N 









E s 


-12.17 


-16.17 






E p 


-7.97 


-8.77 







TABLE I: TB parameters in electron volts. The first nearest-neighbor interactions are given for 
an interatomic distance do (A) 
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a) 




V=0 

FIG. 1: Molecular structure of a) azulene, b) tetracene, c) Ci 6 H 33 Q-3CNQ, d) OHAPy-C-DNB. 
The line at the left side of the molecules defines the zero of the potential \ e xt- 
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FIG. 2: Voltage dependence of the dipole moments along the axis z (straight lines: TB; dashed 
lines: LDA). Inset: dipole moments along z and x axes for OHAPy-C-DNB. 
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FIG. 3: Voltage dependence of the energy levels of tetracene (straight lines: TB; dashed lines: 
LDA). 
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FIG. 4: Voltage dependence of the energy levels of Q-3CNQ (straight lines: TB; dashed lines: 
LDA). 
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FIG. 5: Voltage dependence of the energy levels of OHAPy-C-DNB (straight lines: TB; dashed 
lines: LDA). 
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